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We have searched for — > hh decays, where h stands for a charged or neutral kaon, or a charged 
pion. These results are based on a 23.6 fb" 1 data sample collected with the Belle detector on the 
T(5S) resonance at the KEKB asymmetric-energy e + e~ collider, containing 1.25 million B^*^ Bi*^ 
events. We observe the decay B® — > K + K~ and measure its branching fraction, B(B® — > K + K~) = 
[3.8l Q'g (stat) ± 0.5(syst) ± 0.5(/ s )] x 10 -5 . The first error is statistical, the second is systematic, 
and the third error is due to the uncertainty in the B® production fraction in e + e~ — > bb events. No 
significant signals are seen in other decay modes, and we set upper limits at the 90% confidence level: 
B(B° -> K~tv + ) < 2.6 x 1(T 5 , B{B° -> ty+ty~) < 1.2 x 1CT 5 and B(B° -> K°K°) < 6.6 x 1(T 5 . 

PACS numbers: 13.25.Hw, 14.40.Nd 



The recent observation of a significant difference be- 
tween direct CP violation in B° — > K^it^ and B ± — > 
K^tt was unexpected and has generated much dis- 
cussion. Possible explanations for this difference include 
a large color-suppressed tree amplitude Q, new physics 
in the electroweak penguin loop jfj|, or both [Bj]. Simi- 
lar measurements of charmless two-body B° s decays may 
provide additional insight into this and other aspects of 
B decays. For instance, a comparison of the CP violat- 
ing asymmetries between the B° and B® may discrim- 
inate among new physics models 0; the angles (f>i(/3) 
and 03(7) of the unitarity triangle may be extracted us- 
ing the time evolution of the decays B° — > 7r + 7r~ and 
B° -> K+K- 0; the branching fractions and CP vi- 
olating asymmetries of these two decays provide infor- 
mation on U -spin symmetry breaking [8j ; and the decay 
B n s — > K~n + can be used to determine (^3(7) @. 

The decay B® — > K + K~ is of particular interest be- 
cause its branching fraction is expected to be large, in 
analogy to that of B° —> K + ir~ , and the final state is 
a CP eigenstate. The time-dependent CP asymmetry 
of this decay is sensitive to the B® — B® mixing phase 
(4> s ) and the width difference of the two B® mass eigen- 
states (AT S ); these two parameters provide a clean probe 
of new physics beyond the Standard Model. CDF and 
D0 have performed a time-dependent CP analysis us- 
ing B® — > J/tpcj) events to measure tfi s and Ar s . The 
results are limited by statistics and no significant devia- 
tions from the SM expectation are observed (lfj| . 

Experimental results to date on charmless B° decay 
have been limited to just a few measurements from 



CDF 11- 13| and Belle jl4j ]. In this paper, we report 
on a search for B® decays to K + K~, K°K°, K~ir + and 
7r + 7r~ based on a (23.6±0.3) fb _1 (Lint) data sample col- 
lected at the T(5S) resonance with the Belle detector 
operated at the KEKB asymmetric-energy (3.6 GeV on 
8.2 GeV) e + e~ collider [15J. In an earlier study, half 
of the center-of-mass (cm.) energy was measured to be 
£* cam = (5433.5 ± 0.5) MeV [ijj. At this energy, the 
total cross section for production of light quark pairs of 
the first two families is around 2.446 nb |l7[ while the 
cross section for bb events is c^ 5 ^ = (0.302 ± 0.014) nb, 
of which a fraction f s ~ (19.5^2 3)% contains B® mesons 
[lfjj ]. Three production modes are kinematically allowed: 
B°B°, B* S B° S and B*B*, where the fraction of B*B* is 
f B . B . = (90.li|;^ ± 0.2)% [20]. The number of B*B* 

pairs is thus computed as N B , B , — L ult x cr^ 55 '' x f s x 
f B , s , = (1.25 ±0.19) x 10 6 . ' 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel-like ar- 
rangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to de- 
tect K\ mesons and to identify muons_(KLM). The de- 
tector is described in detail elsewhere [21] . 

Charged kaons and pions are required to have a dis- 
tance of closest approach to the interaction point (IP) 
of less than 3.0 cm in the beam direction and less than 
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0.3 cm in the transverse plane. Charged kaons and pi- 
ons are identified using dE/dx measurements from the 
CDC, Cherenkov light yields in the ACC, and timing in- 
formation from the TOF. This information is combined 
in a likelihood ratio, TZk/tt = £>k/ (At + £>k), where Lk 
(Cir) is the likelihood that the track is a kaon (pion). 
Charged tracks with TZx/n > 0.6 are treated as kaons, 



and with TZk/tt < 0.6 as pions |22|. Furthermore, charged 



p*2 
J bcam 



p* B 2 , and the energy difference, 



SJ eam , where p* B and E* B are the momen- 



tracks positively identified as electrons or muons [22j are 
rejected. With these selections, the kaon (pion) identi- 
fication efficiency is about 83% (88%), while 12% (8%) 
of kaons (pions) are misidentified as pions (kaons). Neu- 
tral kaons are reconstructed in the Kg — > tt + tt~ decay 
channel and are required to have an invariant mass in the 
range 490 MeV/c 2 < M„ +7r - < 510 MeV/c 2 . The inter- 
section point of the tt + tt~ pair must be displaced from 
the IP [iH. 

B® candidates are selected by combining kaons and 
pions in appropriate pairs and separated from back- 
ground using two variables: the beam-energy-constrained 
mass, M bc = y/E\ 
AE = E 

turn and energy of the reconstructed B® meson in the 
cm. frame, respectively. Figure 1 shows the GEANT- 
based |24| Monte Carlo AE-Mb c distributions for the 
B9 x — > hh candidates from various two-body, three- 
body and four-body T(55) decays generated with a B 
meson decaying into an hh pair. Although only one B 
meson per event is fully reconstructed, we can identify 
the T(55) decay from which it originates based on its 
location in the AE-M^ C plane. Candidates with —0.2 
GeV < AE < 0.2 GeV and 5.35 GeV/c 2 < M bc < 5.45 
GeV/c 2 are selected. Since the dominant source of B° 
mesons is T(55) — > B*B*, we search for B® mesons 
only in this decay channel and define the signal region 
to be -0.1 GeV < AE < 0.0 GeV and 5.40 GeV/c 2 
< M hc < 5.43 GeV/c 2 . 

After applying the M^ C -AE selection, there are 
14528, 30613, 27454, and 444 candidates for the 
K + K~, K~ir + , 7r + 7r~ and K°K° modes, respectively. 
These candidates are predominantly from continuum 
events, i.e., e + e~ — > qq, where q stands for a u, d, s 
or c quark. The event topology difference between qq 
and bb events is exploited by computing a Fisher dis- 
criminant 25J based on a set of modified Fox- Wolfram 
moments [26J. Signal (C s ) and background (C q q) likeli- 
hoods are formed using a Monte Carlo (MC) simulation 
and data outside the signal region, respectively. They 
are combined into a likelihood ratio TZ = C s / (C s + C q q). 
The selection criterion, based on TZ, is determined by 
maximizing S/y/S + B, where S and B are the num- 
ber expected in the signal region of signal or back- 
ground events, respectively. The expected signals are 
determined by assuming the following branching frac- 
tions [27]: B{B® ->■ K+K~) = 2.6_x 10" 5 , B(B° 
K-7T+) = 4.6 x 10~ 6 , B{Bl K°K°) = 1.2 x HT 5 , 




FIG. 1: Monte Carlo distributions of AE-M bc for B)- s) -> hh 
candidates from various T(5S) decay modes with B mesons. 
Events in the circles are fromT(55') -» B° W B° W ; candidates 
in the parallelograms are generated withT(5S) -»■ B° M B° M ; 
three-body B^B^-k and four-body BBtytt events are lo- 
cated at M bc > 5.35 GeV/c 2 and AE < -0.05 GeV. 



B(B° 



') = 1.0 x 10" 



For the B\ 



K+K- 



mode, we apply a looser criterion on TZ if the event 
contains an identified electron (muon) with momentum 
larger than 0.5 (0.8) GeV/c. After the TZ requirement, 
300, 444, 188 and 345 candidates are retained for the 
K + K ~, K~ir + , 7r + 7r~, and K°K° modes, respectively. 

Backgrounds from B meson decays are studied using 
large MC samples, which include T(55) -> B^* ] B { s * ] , 
T(5S) -> B*Bir and T(55) BBirir events. _ The 
contributions from T(55) -> BB, T(55) -> B*B and 
T(55) B*B* are negligible since the hh candidates 
from the corresponding B decays lie outside the required 
M-^c-AE region as shown in Fig. [T] Out of the four 
B® decays, B meson backgrounds only appear in the 
B® — > K~ix + mode. A non- negligible contribution from 
T(5S f ) — > B s B s events is present when one of the 
kaons from B® — > K + K~ is misidentified as a pion 
(cross- feed) . The second B meson background is the 
B° -> R-TT+ events from_three-body T(55) -> B*Bn 
and four-body T(5S f ) — > BBirir decays. With the branch- 
ing fractions of T(55) -> B*Bn and T(55) -> BBirir 
assumed to be 6.8% and 9.2%, respectively (28|, we ex- 
pect to reconstruct about five B° —> K~ir + events, lo- 
cated outside the signal region. These cross-feed and 
B° — > A'~7r + backgrounds are considered when extract- 
ing the B° s — > K~ir + signals. 

We perform an unbinned extended maximum likeli- 
hood fit to Mbc and AE to extract signal yields. The 
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likelihood function is defined as 



_V JV- N 



TABLE I: Contributions to the systematic error (%). 



Nl 



(1) 



i=i 



where N is the total number of events, i runs over the se- 
lected events and j over the signal and background com- 
ponents. Nj is the number of events for component j, 
and Pj is the corresponding probability density function 
(PDF). The continuum PDF is the product of a second- 
order polynomial function for AE and an empirical AR- 
GUS function [29j for Mb c . For each mode, the signal 
PDF is modeled from MC with a Gaussian function for 
Mb c and a double Gaussian for AE. The mean values of 
Mbc and AE are calibrated with B Q S — > D+tt~ decays, 
and the AE width is calibrated with D° — > K + tt~ de- 
cays. For the B° -> R-TT+ mode, the B° K + K~ 
cross-feed and the B Q — > K~ir + background are mod- 
eled by two-dimensional smoothed histogram functions. 
Yields for signal and continuum candidates, and the pa- 
rameters of the continuum PDF, are allowed to float in 
the fit while the parameters for other components are 
fixed. The branching fraction (B) is computed as: 



B 



e x 2N 



(2) 



B*BZ 



where N s is the fitted signal yield and e is the MC effi- 
ciency. 

Two types of systematic uncertainties are considered: 
uncertainties associated with the fit and uncertainties on 
the signal reconstruction efficiency and number of B® 
meson pairs. The fit systematic uncertainties are due 
to the modeling of the signal and continuum PDFs, and 
the statistical uncertainties in the background yields that 
were fixed in the fit. The uncertainties due to the sig- 
nal PDFs are obtained by varying each PDF parameter 
successively by one standard deviation and repeating the 
fit. The systematic uncertainty is the quadratic sum of 
the changes in the signal yield. The uncertainty in mod- 
eling the continuum background is studied by changing 
the AE PDFs from second- to first-order polynomials. 
For the B® — > K~ir + mode, the fit is repeated with the 
B° —> K + K~ cross-feed yield varied by plus or minus 
one standard deviation and the signal yield variations 
are assigned as systematic uncertainties. The systematic 
error that arises from the B° — > K~ir + background is 
obtained by taking the difference of the signal yield with 
and without including the B° — > K~n + PDF in the fit. 

The second type of systematic uncertainty is deter- 
mined as follows. For the TZ requirement, we use the 
decay B° — > D~n + to estimate the discrepancy between 
data and MC. The same event selection except the con- 
tinuum suppression used in Ref. [2(| is applied to re- 
construct B® — > D~ir + candidates, where the D~ me- 
son is identified via the D~ — > <f>7r~ , D~ — > K a K~ and 
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D~ — > K*°K~ decays. When forming the variable TZ, 
the Dj mesons are treated as stable particles to mimic 
the B® — > hh events and the same sets of weighting fac- 
tors used to combine the modified Fox- Wolfram moments 
in the hh analysis are adopted. We compare the reduc- 
tion fractions in the D~tt + data and MC with the TZ 
requirements for the four hh modes to obtain the system- 
atic uncertainty. The data-MC differences with various 
TZ requirements are all less than 2/3<r and we conserva- 
tively assign the quadratic sum of the data-MC difference 
and the statistical uncertainty on the Djir + sample as 
the systematic uncertainty. 

The identification of kaons and pions is calibrated us- 
ing a control sample of D* + — » D q (K^tt + )tt + decays. 
For two-body B® — > hh decays, this systematic uncer- 
tainty is 0.7% per kaon and 0.6% per pion. The K s 
reconstruction efficiency is verified using a sample of 
D + —> K s tt + and D + — s- K~tt + it + decays. We com- 
pare the ratio of the yields of the two decay modes with 
the Monte Carlo expectation, which is obtained by gen- 
erating a large Monte Carlo sample with the proper con- 
tinuum and BB fractions. A systematic error of 4.9% 
per K s meson is obtained by adding, in quadrature, the 
deviation of the data and MC ratios and the uncertain- 
ties of the branching fractions of the two decay modes, 
where the latter is the dominant error. The systematic 
uncertainty due to the track reconstruction efficienc y is 
estimated using partially reconstructed D* events [3CJ 
and is 1% per track. Sources of uncertainty in the num- 
ber of B*B* pairs include Lint, j fs, and 
Systematic uncertainties are summarized in Table HI 

The fit results are shown in Figure [2] and summa- 
rized in Table [Til A significant signal is observed in the 
B® — > K + K~ mode, and the branching fraction is mea- 



sured to be B = [3.8±^(stat)±0.5(syst)±0.5(/ s )] x 10~ 5 



with a significance of 5 
fined by £ = ^21n(£ max /£ ), where C 



a. The signal significance is de- 
(£o) is the 
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AE (GeV) M bc (GeV/c 2 ) 



FIG. 2: Distributions of AE (Mb c ) with fit results superim- 
posed for the K + K~ (a,b), K + tv~ (c,d), n + iv~ (e,f), and 
K°K° (g,h) events in the M bc (AE) signal region. The 
blue solid curves represent the fit results, in which the red 
dot-dashed (grey dashed) curves represent signal (continuum 
background). The green dotted curves in the K~tt + plot rep- 
resent the K + K~ cross-feed. 

likelihood value at its maximum (with zero signal yield) 
obtained after convolving the likelihood function with a 
Gaussian function having width equal to the fitting sys- 
tematic uncertainty. For the other decay modes, the 90% 
upper limit (B go %) is computed as 

[® 90% £(B)dB 

^-i — — = 0.9, 3 

Jo £(B)dB 

with the likelihood function after convolving with a Gaus- 
sian width equal to the total systematic uncertainty. 
In conclusion, we observe — > K + K~ with 

B(B° S ->• K+K-) 
= [3.8to;Sj(stat) ± 0.5(syst) ± 0.5(/ s )] x 10~ 5 . (4) 

Our result is consistent with the Standard Model predic- 
tion H and the CDF measurement ([2.44±0.14±0.46] x 



TABLE II: Summary of the signal yields, significances (E), 
reconstruction efficiencies (e), branching fractions (B) and up- 
per limits (U.L.) at the 90% confidence level. 
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10~ 5 ) [12]. No significant signals are observed in the 
other modes, and we set upper limits at 90% confidence 
level: 

B{B° S -> K-TT+) < 2.6 x 10~ 5 , 

B{B° S ->■ 7T+7T-) < 1.2 x 10~ 5 , 

B{B° S ->■ K°K°) < 6.6 x 10~ 5 . (5) 

The first two limits are consistent with results from 
CDF 13j], although with less sensitivity, and the third 
is a first report: this decay is very challenging to recon- 
struct at a hadron collider. 
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